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Abstract 

The lowest-order scaling consideration of the magnetic state formation in the 
Kondo lattices is performed within the s — f model with inclusion of anisotropy 
for both the f — f coupling and s — f exchange interaction. The Kondo 
renormalizations of the effective transverse and longitudinal s — f coupling 
parameters, spin-wave frequency, gap in the magnon spectrum and ordered 
moment are calculated in the case of both ferro- and antiferromagnets. The 
anisotropy-driven change of the scaling behavior (e.g., critical value of g for 
entering the strong-coupling region and the corresponding critical exponents) 
is investigated numerically for N = 2 and analytically in the large- limit. 
The dependence of the effective Kondo temperature on the bare s — f coupling 
parameter g weakens in the presence of anisotropy. The relative anisotropy 
parameters for both the s — f and f — f coupling are demonstrated to de- 
crease during the renormalization process. The role of next-nearest exchange 
interactions for this effect in the antiferromagnet is discussed. 
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I. INTRODUCTION 



Anomalous rare earth and actinide compounds, including so-called Kondo lattices and 
heavy- fermion systems, are studied extensively by both experimentalists and theorists 0-0 . 
It is now a common point of view that the most interesting peculiarities of electronic and 
magnetic propeties of these systems are due to the interplay of the on-site Kondo effect 
and intersite magnetic interactions. Whereas the one-impurity Kondo problem, being itself 
very difficult and rich, is now studied in detail |]^, the Kondo-lattice problem is still a 
subject of many investigations Usually this problem is studied within the standard 

s — f exchange or Anderson models. On the other hand, strong effects of crystal field and 
anisotropic interactions are expected in anomalous 4/ and 5/-systems (see, e.g., 0). These 
effects can lead to anisotropic terms in the Hamiltonian. It is well known that the change 
of symmetry of the s — f exchange interaction modifies qualitatively the infrared behavior 
in the one-impurity case ||10|,|7|,pJ]| . Thus one could expect that similar effects should take 
place in the lattice case. Therefore a question arises whether anisotropic contributions are 
important also in the problem of competition of the Kondo effect and magnetism. It should 
be noted that this question is relevant not only for magnetic systems, but also for models 
with pseudospin degrees of freedom (e.g., for strongly anharmonic crystals demonstrating 
band Jahn- Teller effect 111 211). 



Theoretical investigations of the Kondo-lattice problem use as a rule methods appropriate 



for calculating low-temperature properties in the strong-couping regime (1/A^-expansion [113 



slave-boson technique). However, these methods are not convenient for the description of the 
transition to the weak-coupling regime (in particular, even derivation of the standard Kondo 
logaruithms is here a non-trivial problem ) . In our previous paper we have proposed an 
alternative approach which starts from the weak-coupling regime and is based on summing 
up leading divergent terms by the renormalization group method. We have investigated the 
formation of magnetic state in the periodic s — f exchange and Coqblin-Schrieffer models 
with the /-subsystem being described by the isotropic Heisenberg Hamiltonian. The aim of 
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the present paper is the investigation of formation of the magnetic Kondo-lattice state for 
various magnetic phases with account of the anisotropy in both the locahzed-spin subsystem 
and s — f exchange interaction. 

In Sect. 2 we discuss the theoretical model and calculate the logarithmic Kondo correc- 
tions to the spin- wave spectrum of anisotropic metallic ferro- and antiferromagnets. In Sect. 3 
we derive the lowest-order scaling equations for the effective transverse and longitudinal s — f 
exchange parameters and renormalized magnon frequencies. In Sect. 4 we present a simple 
analytical solution with magnon spectrum renormalizations being neglected, which is possi- 
ble in the large- limit of the Coqblin-Schrieffer model. In Sect. 5 we discuss results of the 
numerical solution of the full scaling equations for = 2 in the presence of the anisotropy 
in localized-spin system only and investigate new features which occur in comparison with 
the isotropic case. In Sect. 6, influence of the anisotropic s — f coupling is considered. 

II. THEORETICAL MODELS AND KONDO CORRECTIONS TO THE 
SPECTRUM OF SPIN EXCITATIONS 

To treat the Kondo effect in a lattice we use the s — d{f) exchange model 

H = Y: tkcLck. + Hf + H,f = Ho + (1) 
where tk is the band energy. We consider the pure spin s — d{f) exchange Hamiltonian with 

Hf = Y: JcS_qSq + E ^^-q^q - KY.iSD', (2) 
q q i 

Hsf = - X! [-^ll'5'k-k'(CktCk't - cj,|Ck/^) + -^±('S'k_k'CkiCk'T + '5'j^-k''^kTCk'T)] (3) 

kk'a/3 

where Sj and Sq are spin operators and their Fourier transforms, r/ > and i^' > are the pa- 
rameters of the two-site and single-site easy-axis magnetic anisotropy in the /-subsystem, re- 
spectively. Note that our consideration can be formally generalized on the Coqblin-Schrieffer 
model with arbitrary (cf. Ref. |^) or to a more general form of the s — f coupling pa- 
rameter matrix [jlT|. For simplicity, we neglect k-dependence of the s — f parameter which 
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occurs in the degenerate s — f model due to the angle dependence of the coupling (see Ref. 
PI). Of course, in fact the f — f exchange has usually the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) origin and is determined by the same s — f coupling. Thus, generally speak- 
ing, the anisotropy of the s — f coupling and /-subsystem are not independent. However, 
crystal-field effects are known to be more important in formation of magnetic anisotropy in 
rare-earth metals than anisotropic exchange interactions fl^ (in this case, the anisotropic 
s — f coupling is obtained by expansion in the parameter kpTj, Vf being the /-shell radius, 
and contains, unlike (^), terms of the type (kS)(k'S)). On the other hand, the situation, 
where anisotropy occurs in the s — f coupling only, may be also considered: this corresponds 
to the strong "direct" / — / exchange (e.g., superexchange) interaction which is character- 
istic for some /-compounds. In the Coqblin-Schrieffer model, which is more appropriate for 



cerium compounds, crystal field results in occurrence of anisotropic s — f coupling |1T5|JT6 
and new excitation branches 1^,^. For simplicity, we restrict ourselves to treatment of a 
single magnon mode in the simplest s — f model (|1|). 

In the ferromagnetic (FM) state the spin-wave spectrum for the Hamiltonian (^ reads 

UJq = UJo + Uexiq), (4) 

= 25(Jq - Jo),^o = -2S7]Jo + (25 - 1)K (5) 

To find the Kondo logarithmic corrections to the spin-wave spectrum we calculate the 
magnon Green's functions in the model (|l|). For a ferromagnet we obtain to second or- 
der in I (cf. the calculations in the isotropic case [|18| ) 



((^ql^i)).. = [uJ-uj^-2 ^( Jo + Jq_p - Jp - Jq + a;o/25) (6j,6p) 

- -1 

(6) 



p 

1 



2. ) 



where we have taken into account kinematic requirements in the magnon anharmonicity 
terms by introducing the factor of (25 — 1)/2S at K (this replacement may be justified 
by considering higher-order terms in the formal parameter 1/S), nk = n{ti^) is the Fermi 
function, 



(7) 



to + tk — tk+p-q — 

The averages that enter can be obtained from the spectral representation for the Green's 
function (|^) to first order in l/S* and contain the singular contributions 



(8) 



with 



^FM ^ ^k(l - »-k+q) 

^(tk-tk+q-^,^2 



(9) 



1^. "K+q "^q; 

Expanding the denominators of in the magnon frequencies we obtain the singular cor- 
rection to the pole of the magnon Green's function 



-2S{ll + Jj) Jp - Jq_p + Jq - Jo + UJo/2S)<^ 



FM 
P 



This result can be represented as 



sumMcl) = -ill + )(i - «q) E K 
W^o = -(/! + /f)E'^J 



FM 



^FM 



(10) 

(11) 



where < aq < 1. Passing into real space (see p[) yields 



E^R (Ot.=i^. [l-cosqR]/^JK[l-cosqR], 



(12) 



R R 

In the approximation of nearest neighbors at the distance d, the quantity a does not depend 
on q. For a spherical Fermi surface we have 
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t^—Ec 



sin kpd 



(13) 



Hereafter we put a = const. Then we may use in further consideration of the scaling 
equations a single renormalization parameter, rather than the whole function of q. 

Now we consider a two-sublattice antiferromagnet (AFM) with the wavevector of mag- 
netic structure Q, 



{S!) = ScosQR„ {S!)= {3^ = 

{Jq = Jmin < 0; 2Q is equal to a reciprocal lattice vector, so that cos^ QRj = 1, sin^ QR^ = 
0; only in this case we can retain the definitions of I± and I\\ in the local coordinate system). 
Passing to the spin-deviation operators in the local coordinate system where 

= St cos QR„ = Sf cos QR„ = (14) 

we derive 

Hf = const + EtCq&i&q + ^D^ib-nK + K^-n)] + - (15) 

Cq = 5( jQ+q + Jq - 2 Jq(1 + 7])) + {2S - 1) K , = S{J^ - Jq+q) (16) 

Diagonalizing ( pTSf ) we obtain the spin-wave spectrum 

^q = Cq-^q^^O+^L(q) (17) 
U;,M = 25(Jq - JQ)^/'(JQ+q " JqY/^ (18) 

= 25( Jo - Jq) [{28 - \)K - 2SriJQ] (19) 

where we have neglected a weak wavevector dependence of uq. 

The Kondo correction to the spectrum ([T7|) reads (cf. Refs. []8|,p!9[|) 

6u;l = -2j2[llujl + 2S^ll{Jci+^ + Jq - 2Jq)( Jp + Jq+p - Jq+q^p - Jq_p) (20) 

+2(/J - /i)(CqCp_q - /^q/^p-q)]$^™ 

with 

^AFM ^ 77.1^(1 — TT-k+q) , ^ 

k l^k - Ik+qj - t^q 

For an antiferromagnet in the nearest-neighbor approximation (Jq+q = — Jq) we obtain 
from (|20|) 

ScolM/ ^el(q) = -2[ll - «(/|f - /!)] E (22) 



^CoVc^o^ = -2[/i - (1 - a).(/f - ID] E ^r^" 



p 



where s = AS'^J^/u^ > 1. 

Introducing next-nearest-neighbor interactions and putting for simphcity /y = 1^ = 1 
we obtain 

5ul = -2[col - 85^«f (Jf - - 4^)] E <^ (23) 

p 

where 

correspond to the contribution of nearest and next-nearest neighbors, and a^^ is given 
(0) with the sum over the next-nearest neighbors. Provided that next-nearest neighbor 
exchange interaction is ferromagnetic (so that the AFM ordering is stable), J^'^ — Jq^ > 
and the next-nearest neighbors result in decreasing the Kondo suppression of the magnon 
frequency, as well as nearest neighbors in the FM case. Then, instead of (p2|), we can use 
phenomenologically (e.g., in the long- wave limit) the expression 

Su^eM/ ^eM = -(1 - E *P™ (24) 

P 

with a' (X a 

(2)j(2)/j(i)_ opposite case of AFM next-nearest neighbor exchange the 

situation is more complicated. In particular, the simple collinear antiferromagnetism can 
become unstable, and formation of the spiral structure is possible. Note the difference 
between the FM and AFM cases by a factor of 2, which is due to violation of time-reversal 
symmetry in a ferromagnet (terms that are linear in uj give a contribution). 

The quantities (^, (|2l|) determine also the singular correction to the (sublattice) mag- 
netization 

6S/S = E = -II E ^r'^^"" (25) 
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III. SCALING EQUATIONS 



Using the perturbation theory results we can write down the system of scahng equations 
in the case of the Kondo lattice for various magnetic phases. Their derivation in the isotropic 
case is described in detail in Ref . . As well as in this paper, we apply the "poor man scal- 
ing" approach ||2D|. In this method one considers the dependence of effective (renormalized) 
model parameters on the cutoff parameter C {—D < C < 0, here and hereafter the energy is 
calculated from the Fermi energy Ep = 0) which occurs at picking out the Kondo singular 
terms. 

The renormalization of /|| is obtained from renormalization of the magnetic splitting in 
electron spectrum, and of Ij_ from renormalization of the second-order contribution to the 
electron self-energy (see corresponding perturbation expressions for a ferromagnet in Ref. 



21|). The renormalized /_|_ chosen in such a way coincides with the three- leg vertex (with 
two electron lines with opposite spins and one magnon line) which yields the most natural 
definition in a magnetically ordered case and agrees with the one-impurity scaling consider- 
ation To find the equation for the effective coupling parameters I^fiC) {I"j{—D) = 1"^) 
we have to calculate the contribution of intermediate electron states near the Fermi level 
with C < tk+q < C + 6C in the sums that enter expressions for the self-energies (which 
include, unlike Ref. [§, magnon frequencies with a gap). Then we obtain 

Sl!r{C) = 2pllv{-^,-^)6C/C (26) 



5Ji,(C) = 2p/^/||r^(-^,-^)5C/C 



where p is the density of states at the Fermi level, uj^x is a characteristic spin-fluctuation 
energy, lo^ is the gap in the spin- wave spectrum, //(x) is a scaling function which satisfies 
the condition 7^(0) = 1 which guarantees the correct one-impurity limit. For both FM and 
AFM phases we have 

FM,AFM( ^ex ^On_//, 2 /r^2N-l\ /QyN 
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where the magnon frequencies are given by (|T^. The corresponding analytical expres- 
sions are presented in Appendix. 

The C-dependent renormalizations of the spin-wave frequencies and ground-state mo- 
ment are obtained in the same way as in the isotropic case from (|T0|), pO| ) (^2]), 
(P5|) and expressed in terms of the same scaling functions. Introducing the dimensionless 
coupling constants 

(we will drop sometimes the index ||, but not ±, so that g^fiC) = glf{C)) and replacing 
Qa 9efi^)^ ^ex{C), ^0 ^o{C) in the right-hand parts of (^) and expressions for 

6Uex{C),Sujo{C) and 6Sef{C), we obtain the system of scaling equation 

dgl{C)/dC = -[g^jiC)rA (28) 
dg^f{C)/dC = -glf{C)g^f{C)A (29) 

d\nuJex{C)/dC = aA/2 X i (30) 

[[g^fiC)? - a{[glf{C)f - [g^j{C)r} AFM 

.ln.o(C.)/aC^.A/2x(^[^-(^)]^^[^-(^)]^>/^ ™ (31) 

[[g^fiC)]' - sil - a){[glf{C)r - [g^fiC)]} AFM 

dlnS,f{C)/dC=[g^f{C)]'A/2 (32) 
where 

A = A{C,u;UC),uo{C)) = r]i~uJ,fiC) /C, -uo{C)/C)/C, (33) 

and 

2(1 -a) FM I 2 FM 

a={ ' ,b={ (34) 

1 - a' AFM 1 AFM 



The integral of motion of the system (^81) , (^91) reads 

[glfiC)]' - [g^fiC)r = f^' = gl-gl = const (35) 



so that the equation (pS]) takes the form 

^g,AC)/^C = -[g%{C)-^,']^ 



(36) 



IV. ANALYTICAL SOLUTION IN THE LARGE-iV LIMIT 

In the formal large- limit in the Coqblin-Schrieffer model where 2 — > in Eqs.(|30D- 
(^)) we can neglect renormalizations of magnon frequencies (note that the true large- iV 
limit in the FM case is somewhat different since symmetry of spin-up and spin-down state 
contributions is violated for N > 2, see Ref. Q). Note that the same approximation is valid 
for N = 2 provided that g is well below the critical value gc for entering the strong-coupling 
region. Then, on taking into account (|A1| ), (|A4| ), equation (|36D can be integrated analytically 



to obtain 



-[arctanh(^/^e/(C)) - arctanh(/i/^)] = G{C) = - ~C^^^~^' ^^^^ 
^^(C) = \n\C/D\ - ((l + w)/2)[(C/cJ- l)ln|l -cJ/C| - {C/uJ + l)ln\l + uJ/C\] (38) 



+ {w/2)[{C/u;o - 1) In |1 - iJo/C\ - (C/c^o + 1) In |1 + cuo/C|] - 1 

G^fiHC) = In \C/D\ - 1[(1 + w')iCyuj' - 1) In 1 1 - u'/C'l (39) 

-w^C^uj^o - 1) In |1 - ujI/C^\ + 1] 
,1 



GlT(C) = [9{C' - uj') + 9{u;', - C')] H-{^\C^ - + ^\C^ - u;'\)/D) (40) 
where 



UJ 



uol + AFM 



(41) 



The scaling trajectories described by (p7|)-(|40|) are shown in Fig.l for yU 7^ 0, = and in 
Fig. 2 for w 7^ 0,/i = 0. Note that these pictures describe adequately the case N = 2, since 
g = 0.15 is considerably lower than the critical values. 
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The anisotropy of s — / coupling results in that the dependence l/(7e/(0 becomes non- 
linear at small ^ = In |D/C| where the one-impurity behavior takes place 

- /xtanh[arctanh(/i/^) - fiQ (42) 

However, this non-linearity is not too strong even for fi/g = 2/3 (Fig.l). Of course, the 
curves with /i 7^ go considerably higher, since the bare coupling parameter g± decreases 
with /I. 

For w the function 1/ gef{C) has a minimum both in the AFM and FM cases, position 
of which, Cinin, is determined by ([A^)- (|X7|) (in the isotropic case, the minimum occurs in the 



3D AFM case only). The minimum may result in non-monotonic temperature dependences 
of physical quantities which are sensitive to the Kondo effect, e.g., of the effective magnetic 
moment. The depth of the minimum A = Gmin — is given by 



A 



\[{l + w)\n{l + w)-w\nw] FM 
i [In 2 + (1 + ti;2) ln(l + w'^) - In w^] 3D AFM (43) 



\n{w + VI + ^2) 2D AFM 

Note that at if ^ 1 we have in all the cases A oc Inw. 
The critical value Qc is determined by the condition 

l/g^ = -(l/fi) tanh(/iG'min) - |G^in - (fi'^ /3)GiJ, fi <^ 1 (44) 

where 



— Gr^ 



\ + l + l[wlnw + {l-w)\n{l + w)] FM 

A + |[1 + ln2 + ln(l + «;2)] 3D AFM (45) 
A 2D AFM 



In the FM case anisotropy in /-system results in an increase of Qc (|G'(0)| increases with w 
more rapidly than IGminl), in the 3D AFM case Qc decreases, and in 2D AFM anisotropy does 
not influence Qc- The effective coupling constant g* = gef{0) remains finite aX g —* Qc — 
and tends to 
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g* = fi/ tanh(/iA). 



(46) 



The effective (renormalized by spin dynamics) Kondo temperature is determined by the 
condition 1/ gef{—T^) = and satisfies the equation 

- G{-Tl^) = arctanh(/i/^) (47) 

Due to the minimum, is also finite at g ^ Qc + 0, = |Cmin|. 

V. EFFECTS OF THE F-SYSTEM ANISOTROPY ON SCALING BEHAVIOR 



Now we treat the physically real case N = 2 with the anisotropy being present in the 
/-system only. To this end we have to consider the full scaling equations for fi = 0,w 0. 
The most important circumstance to be taken into account is the renormalization of the 
magnon frequencies. Writing down the first integrals of the system (pHD, ( pO] ) 



QefiC) + (2/a) InuJexiC) = const 
gefiC) + (2/6) \nuo{C) = const 



(48) 



results in 



uJex{C) = uJexew{-a[9ef{C) - g]/2) 
uJoiC) = uo exp{-b[gef{C) - g]/2) 



(49) 



As follows from (H) 



1/a 



Up 



S 



Substituting (|^) into (^8|) we obtain 



(50) 



d{l/gef)/d^ = -T/(exp(^ - \ex - a[gef - ^]/2), w;exp(^ - \ex - b[gef - g]/2)) (51) 



where 



^ = ln\D/C\, Xex = lia{D/uJex) > l,w = Wo/cUe 



12 



For an antiferromagnet in the nearest-neighbor approximation we have a = b and the equa- 
tion (^) takes the form 

d{l/gef)/d^ = -vl/(Ae. - e + a[gef - g\l2)) (52) 

For finite values of N the singularities of the scaling functions, that occur in the mag- 
netically ordered phases, play the crucial role due to peculiar properies of the differential 
equation (|5l|) . In particular, one can prove p that g^f diverges at some ^ at arbitrarily small 
g (i.e. gc = 0) unless the singularity cutoff is introduced. To make the value of gc finite one 
has to cut the singularities. This may be performed by introducing small imaginary parts, 
i.e. by replacing in (|Al|) , {\K^) 

ln|l-a;| ^ Re ln[l - x(l + = ^ ln[(l - x)^ + (53) 

(1 - x)-'/^9{l - a;) {[(1 - xy + S^x^f^ + 1 - x]/2}VV[(l _ xf + d^x']'/^ 

The x-dependence of the cutoff parameter can be in principle neglected, as in Ref. [Q, 
since it does not influence appreciably numerical results (since 6 is important at x = 1 
only); however, this dependence is needed to pass correctly to the limit ujq 0). As 
one can see from (|2^) , the value of 6 should be determined by the magnon damping at 
q = |k — k'l ~ 2k p. This damping is due to both exchange and relativistic interactions. 
Hereafter we put in numerical calculations 6 = 1/100. We accept also A = \n{D/uj) = 5. 

The dependences gefiO ^r 3D FM and AFM phases according to (|5ll) , (^) at g close 
to the critical value are shown in Fig.3. As well as for the isotropic case |^, there occur 
large intervals of a non-Ferm-liquid (NFL) behavior where magnon spectrum becomes soft. 
In this regime, the relation between the arguments of the scaling function ri{x,y) in ( |5lD is 
flxed by the singularity point \C\ = uj{C) (see Appendix, a;(^) is deflned in the same way as 
in (1^)) or by the condition 

ln[^/^(0] = e - A (54) 
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After substituting ( ^9]) into ( ^41) we see that the dependence gefiO is hnear in ^ only in the 
case a = b where 



r?e/(0-^7-2(e-A)/a 



(55) 



Thus for the FM phase the dependence QefiO NFL region is different from the isotropic 

case. 

In the 2D AFM case the anisotropy does not practically influence the behavior gef{C) 
because of strong singularity of the scaling function. 

The dependences l/g*{g) and ^*{g) ( C,* = ln{D/T^)) for a 3D anisotropic ferromagnet 
are shown in Fig. 4. These dependences are more similar to those in the isotropic antiferro- 
magnet rather than isotropic ferromagnet (see Fig. 5 of Ref. In particular, a wide plateau 
with ^*{g) — ^* can be seen in Fig. 4, whereas in the isotropic ferromagnet ^*{g — > (^c + 0) 
increases in a not too narrow region. The difference is connected with the absence of the 
scaling function maximum in the latter case. Thus the anisotropy makes the dependence 
the effective Kondo temperature on the bare coupling parameter still weaker in comparison 
with the isotropic case. Due to the minimum of the function l/g^fi^), g* remains finite at 
g ^ gc — 0, as well as in the limit N^oo, g*^g* = l/A with A given by (^5]). It should 
be noted that the functionl/5f*(5f) (Fig.4) approaches A at very small \g — gd which are 
practically unreachable. 

On the other hand, strictly speaking, C,*{g) diverges at g ^ gc + 0. The increase of C,*{g) 
takes place also in an extremely narrow region (of order of 10~'^-10^'^) only and is not shown 
in Fig.4. As demonstrate numerical calculations, in this region we have 



The "critical exponents" 7 turn out to be non-universal, decreasing with increasing w, i.e. 
the minimum depth; for the 3D ferro- and antiferromagnets with w = 0.3 we have 7 ~ 0.2 
and 7 ~ 0.05 (the corresponding values in the isotropic case are 7 = 1/2 and 7 ~ 0.1 
respectively 



eig)-C 



iMia - gc)/g],T;^ ~ (g-gc) 



7 



(56) 
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A comparison of the critical parameter values in the isotropic and anisotropic cases is 
presented in the Table 1. One can see that ior N — 2 the anisotropy results in a decrease 
of the critical value gc in all the cases, unlike the large-A'" limit. This decrease is more 
appreciable in the FM phase where linear terms in the anisotropy parameter enter the 
equations. 

Table 1. The critical values gc and ^* for different magnetic phases in the isotropic and 
anisotropic cases &t N — oo (from analytical results, see Sect.4) and N — 2 (from numerical 
solution). The parameter values are A = 5, a = 1/2, a' = 0. For N — 2 , the cutoff 5 = 1/100 
is used and the "critical value" ^* is estimated from the plateau in the dependence Ci.9) 
(see the discussion in the text). 





w 




FM 


3D AFM 


2D AFM 


^ oo 





9c 


0.167 


0.171 


0.176 




0.3 




0.169 


0.170 


0.176 







e: 


oo 


5.35 


5 




0.3 




5.73 


5.31 


5 


N ^2 







0.139 


0.132 


0.127 




0.3 




0.133 


0.130 


0.127 







c 


6.13 


6.07 


6.07 




0.3 




6.24 


6.04 


6.06 



The value of ^* decreases with anisotropy in the AFM case but increases in the FM case. 
This is explained by that in FM cuq is renormalized stronger than uJex, so that co{^) tends to 
zero more rapidly. As demonstrate numerical calculations, ^* decreases with w in the AFM 
case too for sufficiently large a'. 

The experimentally observable quantities can be obtained by using the formulas 

T*=L>exp(-r) (57) 

{g > gc) and 

S* = ScfiC = 0) = SeM-[9* - ^]/2) 
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ul^ = UexiC = 0) = ujexexp{-a[g* - g]/2) (58) 
Uq = uJo{C = 0) = uJoexp{-b[g* - g]/2) 

{g < gc). In particular, we obtain the relation 

f (S*/Sf FM 

^oV^o = (59) 
[ S*/S AFM 

The quantities ( |58D are finite at g gc — 0. However, as follows from the standard scaling 
treatment of the phase transition, in some region we have the law 

uj* = uj(C = 0)r^igc-gr (60) 

Renormalization of relative anisotropy parameter is given by 

^(C) = = eM-^-^[9ef{C) - g]) (61) 

This is shown in Fig. 5. The corresponding temperature dependence can be obtained by the 
replacement \C\ ^ T. 

VI. EFFECTS OF ANISOTROPIC S- F COUPLING 

Strictly speaking, in the case where n ^ the full system of scaling equations cannot 
be simplified. However, a simple analytical treatment is possible in the case ^ g which is 
physically real for magnetic systems. Under this condition we can expand 

aif = ^/gfr-^ - 9ef - ^f^Vgej (62) 

Provided that the expansion holds at ^ = 0, this will hold with increasing g^j too. Then we 
obtain from (p^) , (pO|) , (^T]) the integrals of motion 

1/2 FM 

(7e/(C) - r/iV(7e/(C) + (2/a) lnaJ,,(C) - const, r = <| (63) 

-a AFM 

1/2 FM 

gefiC) - Qii^lgeAC) + (2/6) lna;o(C) ~ const, B = { (64) 

-(1 - a)s AFM 
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so that 

u,^{C) = uJe.^M-ci[9ef{C) -g- T^l\l/g,f{C) - llg)]l2) (65) 
uJo{C) = o7e.exp(-6[(7e/(C) -g- eii\l/g,j{C) - l/g)]/2) (66) 

On substituting (|65| ) , (|66|) into ( ^6l) we obtain the closed equation 

dgef/d^ = -{glf - /i)r/(exp(^ - A^^ - a[gef -g- T^'^{l/g^} - l/^)]/2), 

weMi - - b[gef -g- 0l^\l/9ef - l/^)]/2)) (67) 



Presence of the terms, that are proportional to /i, in the scaling function arguments in 
results in a weak smearing of the linear dependence gefiO NFL regime even for a = b. 

In the FM case, to accuracy accepted, the expressions (^),(^) can be represented as 

oJex(C) = uJe,exp{-a[g^jiC) - gi_]/2) (68) 
uoiC) = exp(-6[4(C) - g,_]/2) (69) 

Further, as follows from (0), in all the cases 

SefiC) = S exp{-[g,j{C) - g]/2) (70) 



Thus the first of relations (^Dj) is violated in AFM case, and the second relation in both FM 
and AFM cases. 

However, the violation owing the anisotropic s — f coupling is weak: renormalization of 
Uq for AFM, which is most appreciable, yields in the exponent of (^) the quantity of the 
order of fi'^s/g oc g only. Thus the most important effect of /i is the deformation of scaling 
trajectories at not too large values of ^, which was considered in Sect. 4 (Fig.l), and main 
corrections to gc are described by 



The case of not small (which can be relevant for pseudospin systems |T2|) can be 
investigated by numerical solving the equations (|28|)- (|3TD. The results are shown in Fig.6. 
One can see that the anisotropy of s — / coupling leads to an increase of the gc values in 
comparison with the isotropic case (cf. Table 1). This is due to both decrease of g± and 
more weak renormalization of uJex according to (|30|) . 
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VII. CONCLUSIONS 



In the present paper we have generahzed the scahng treatment of Ref. ||^ by including 
the anisotropy in both the s — f couphng and /-system itself. We have demonstrated that 
the magnetic anisotropy modifies considerably the scaling behavior in the Kondo lattice 
problem. 

In all the cases, the anisotropy in the /-subsystem (the gap in the magnon spectrum) 
results in occurrence of a non-monotonous dependence (of a maximum) of the effective 
coupling constant QefiCj- This prevents the increase of C{9) at — > (?c + 0, which becomes 
practically not observable even for a ferromagnet (Fig.4), unlike the isotropic case. The 
dependence of the effective Kondo temperature on the bare s — f coupling parameter becomes 
weaker in the presence of the anisotropy. Further, the minimum in the dependence 1/ gef{C) 
results in that g* (and, consequently, cU*^, tUg, and S*) are always finite aX g ^ Qc — Q- As for 
quantitative changes, the anisotropy favors a non-magnetic Kondo state (the critical values 
of bare coupling constant decrease) The critical region of magnetic instability becomes more 
narrow (especially in the FM case), so that the non- Fermi-liquid behavior is suppressed. 
At the same time, anisotropic s — f coupling influences noticeably the form of the scaling 
trajectory for small ^, but becomes not important with increasing gefiO- 

Owing to a more rapid Kondo renormalization of the gap, the system tends to "isotropic" 
behavior of gefiO approaching the strong coupling regime. Such a renormalization (|6lD 
may be important for analysis of the spin excitation spectrum in anomalous /-systems. One 
can expect that the observable renormalized spectrum gap in these systems is relatively small 
and strongly temperature-dependent in comparison with the "usual" /-electron magnets. As 
demonstrate our calculations, next-nearest-neighbor exchange interactions in AFM phase are 
important for this effect. The dependence uJo{T) can be investigated not only by the neutron 
scattering, but also by simple methods like the ferromagnetic resonance. 

The change of the critical exponents of the phase transition at g g^ with changing 
the bare anisotropy parameter w turns out to be very strong. It is interesting that their 
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values depend continuously on the anisotropy and are non-universal (possibly, higher order 
contributions to the scaling equations will change this result; this question needs further 
investigations). It should be noted that, unlike the one-mpurity Kondo problem (where the 
"phase transition", connected with disappearance of the local moment, exists for N = oo 
only, and a crossover takes place for finite N, see Ref. [0]), the phase transition in the Kondo 
lattices is physically real, being a magnetic-nonmagnetic transition. The situation is similar 
to the onset of magnetism in itinerant electron systems. Of course, one has to bear in mind 
that the treatment of this transition within the lowest-order scaling may be only qualitatively 
correct, and a more detailed (e.g., numerical renormalization group) consideration is needed 
to describe the crossover region. 

Similar to Ref. [§, at approaching the critical value of magnetic instability Qc, the tran- 
sition to an "incoherent" regime, where non-spin-wave excitations of /-system play the 
dominant role, should be considered. In this regime, the minimum of l/(7e/(0 can be sup- 
pressed. However, the use of the model "paramagnetic" scaling function for describing the 
incoherent contribution (as in Ref. |Q]) seems to be unreasonable in the presence of the 
anisotropy, since spin dynamics in the paramagnetic state of /-systems in a strong crystal 
field is rather complicated [P,P^. 

The work was supported in part by Grant 96-02-16000 from the Russian Basic Research 
Foundation. 



APPENDIX A: APPENDIX. SCALING FUNCTIONS IN ANISOTROPIC 

MAGNETS 



Using the long- wave approximations uj^J^(q) ~ g^, ujfj^^(q) ~ g in the whole Brillouin 
zone (which is justified, e.g., at small kp), we get from ( pTj) 



x,y) = — In 



1+x + yl— y 



1 — X — yl + y 



(Al) 



where ujex = ^exi'^kp)- For an antiferromagnet we derive 
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1 - r V 1 - 

where ri^^^^{x) is the corresponding scahng function in the isotropic case, 

AFM/ N _ / -X ^\n\l - x^\, d = 3 , 



^(x) being the step function. Then we have 



4 In 



±zlL 



2 



d = 3 



(A4) 

One can see that the logarithmic singularities of the functions 7]^^^{x) = ri^'^^{x,0) and 
?7^^^^(x) at X = 1 are shifted to the points x + y = 1 (FM) and x'^ + y'^ = 1 (AFM), i.e. 
\C\ = ZJ with u defined by (^T]). Besides that, the anisotropy results in occurrence of the 
second singularity at y = 1 (i.e., \C\ = Uq). 

Presence of such singularities is a general property which does not depend on the spec- 
trum model p. The functions ri^'^^{x,y) and r)^^^^ {x^y) {d = 3) change their sign at 
y(x + = 1, i.e. 



\C\ = |Cmin| = UJo{uJo + UJex) (A5) 



and + 2y^ = 2, i.e. 



^7min| = V^0'+^L/2 (A6) 



respectively. For d = 2 the function rj^^'^'^^x) has strong square-root singularities at x^+y^ 
1 and y = 1, and vanishes in the interval 1 — x^ < < 1, i.e. 



ICininI = ^ui + uJl^ > \C\ > Uq, (A7) 

changing its sign after passing this interval (but a smooth contribution occurs for more 
realistic models of magnon spectrum). 

Note that in the limit of strong magnetic anisotropy uJo/uJex oo the singularity at 
y — > 1 becomes very strong: 
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r^\y) = {I - y^)-^ (A8) 

It should be noted that the influence of magnetic anisotropy on spin dynamics can be also 
considered for the paramagnetic phase (i.e., for the problem of the local moment formation). 
In this case the singularity ( [A8[ ) also becomes weaker. If we accept, as in Ref. |p, the "spin 
diffusion" approximation with the spin spectral density 

1 Vq"^ 

-K [uo - uoqY + {Dq'^)^ 

{T> is the spin diffusion constant, uJex = 4T>kp) we obtain 

r)^^(x,y) = — farctan h arctan — - — | (AlO) 

2a; V 1-2/ 1 + 

This function has a flnite jump at \C\ = ujq. Of course, this approximation can be hardly 

justifled for real /-systems. Therefore we do not present concrete calculation results for the 

PM case. However, one can expect that qualitative effects of anisotropy are similar to those 

in the magnetically ordered phases. 

Figure captions 

Fig.l. The dependence l/gej on ^ = \n\D/C\ in the case of anisotropic s — f coupling 
for a 3D antiferromagnet (upper solid line), 2D antiferromagnet (short-dashed line) and 3D 
ferromagnet (long-dashed line). The parameters are A = ln{D/uj) = 5, g = 0.15, /i = 0.1. 
The lower solid line shows the curve for the 3D antiferromagnet with /i = 0. 

Fig. 2. The dependence 1/gef on ^ = In \D/C\ in the presence of anisotropy in /-system 
for a 3D antiferromagnet (solid line), 2D antiferromagnet (short-dashed line) and 3D ferro- 
magnet (long-dashed line). The parameters are A = 5, g = 0.15, w = 0.3. 

Fig. 3. The scaling trajectories gefiO in an anisotropic ferromagnet for g = 0.1333 > gc 
(upper solid line) and g = 0.1332 < gc (lower solid line), and 3D antiferromagnet for 
g = 0.1302 > gc (upper dashed line) and g = 0.1301 < gc (lower dashed line), w = 0.3. 

Fig. 4. The dependences l/g*{g) for g < gc and ^*{g) — A for g > gc in an anisotropic 
ferromagnet with A = 5, a = 1/2, w = 0.3, 6 = 1/100. The dashed line is the curve 1/g — X. 



21 



Fig. 5. The scaling behavior of the effective anisotropy parameter tye/(0 — ^oiO /^ex{0 
for an anisotropic ferromagnet with the same parameters as in Fig. 2 (sohd hnes) and anti- 
ferromagnet with a' = 0.4, g = 0.1375 > Qc (upper dashed hne) and g = 0.1374 < gc (lower 
dashed hne). 

Fig. 6. The dependences l/(7e/(0 ^^e case of anisotropic s — f coupling (/x = 0.1, w — 0) 
for a 3D ferromagnet with g = 0.1619 > gc (lower solid line) and g = 0.1618 < gc (upper solid 
line) and a 3D antiferromagnet with g — 0.1563 > gc (lower dashed line) and g — 0.1562 < gc 
(upper dashed line). 
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